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Background: Well-validated models of maternal behavior in small-brain mammals posit a 
central role of oxytocin in parenting, by reducing stress and enhancing the reward value 
of social interactions with offspring. In contrast, human studies are only beginning to gain 
insights into how oxytocin modulates maternal behavior and affiliation. 

Methods: To explore associations between oxytocin receptor genes and maternal 
parenting behavior in humans, we conducted a genetic imaging study of women selected 
to exhibit a wide range of observed parenting when their children were 4-6 years old. 

Results: In response to child stimuli during functional magnetic resonance imaging 
(fMRI), hemodynamic responses in brain regions that mediate affect, reward, and social 
behavior were significantly correlated with observed positive parenting. Furthermore, 
single nucleotide polymorphisms (SNPs) (rs53576 and rs1042778) in the gene encoding 
the oxytocin receptor were significantly associated with both positive parenting and 
hemodynamic responses to child stimuli in orbitofrontal cortex (OFC), anterior cingulate 
cortex (ACC), and hippocampus. 

Conclusions: These findings contribute to the emerging literature on the role of oxytocin 
in human social behavior and support the feasibility of tracing biological pathways from 
genes to neural regions to positive maternal parenting behaviors in humans using genetic 
imaging methods. 
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INTRODUCTION 

Dysfunctional maternal parenting during early childhood is a 
robust risk factor for mental and physical disorders associated 
with mortality in humans (Wegman and Stetler, 2009). In spite of 
its public health importance, however, the neurobiology of dys- 
functional maternal behavior in humans is poorly understood. 
Research with rodents highlights the involvement of oxytocin 
(OT) in the initiation and expression of maternal approach 
behavior, affiliation, and attachment with young (Lim and Young, 
2006). Such behaviors correlate with both peripheral measures of 
OT in plasma and the expression of oxytocin receptors (OXTR) 
in the brain (Francis et al., 2002). In line with these observa- 
tions, the OXTR gene single nucleotide polymorphism (SNP), 



rs53576, in humans, has been associated with a broad range of 
social behaviors potentially relevant to adaptive parenting behav- 
ior (Costa et al, 2006; Israel et al, 2009; Rodrigues et al, 2009). 
Further elucidating the associations between oxytocinergic poly- 
morphisms and both behavioral and neural phenotypes may 
provide complementary evidence for the involvement of OT in 
human parenting. Here we combine functional magnetic reso- 
nance imaging (fMRI) and genetic analyses to trace biological 
pathways from genes to neural circuits to human maternal par- 
enting. We examine whether the same OXTR polymorphisms 
are associated both with variations in observed maternal par- 
enting and variations in maternal neural responses to child 
stimuli. 
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A consensus regarding neurobiological mechanisms of mater- 
nal behavior in small-brain mammals has emerged (Brunton and 
Russell, 2008; Numan and Woodside, 2010). In such animals, 
being a mature female does not provide a sufficient basis for 
maternal behavior. Rather, nulliparous females often avoid and 
even attack newborns (Numan and Insel, 2003). Current mod- 
els posit that rising levels of estrogens during pregnancy "prime" 
key brain regions, initiating a process that reverses the valence of 
infant stimuli from aversive to rewarding, largely by upregulat- 
ing OT receptors (Carter, 1998; Numan and Woodside, 2010). 
Much of the knowledge regarding the ability of OT to regulate 
social interactions is based on data from animals using cen- 
trally administered agonists and antagonists or knockout mice. 
OXTR gene knockout mice are profoundly deficient in mater- 
nal behavior (Ragnauth et al., 2005; Takayanagi et al., 2005). In 
rats, intra-cranial injection of OT following estrogen priming 
induces maternal behavior in virgin females (Campbell, 2008). 
Reciprocally, an OXTR antagonist inhibits the natural postpartum 
onset of maternal behavior (Insel, 2010). Although initial studies 
suggest similar social effects on humans (see Bartz and Hollander, 
2006, for a review), the relevance of these animal findings for 
humans is not yet clear. 

Oxytocin can function as a synaptic neurotransmitter by 
being released diffusely into extracellular fluid from dendrites, 
influencing distant neural sites (Landgraf and Neumann, 2004; 
Ludwig and Leng, 2006; Ross and Young, 2009; van den Burg 
and Neumann, 2011). Recent morphological and functional data 
also demonstrate that OT neurons emanating from the hypotha- 
lamus have long-range axonal projections capable of releasing 
OT (Ross et al., 2009; Knobloch et al, 2012). OT directly and 
indirectly influences sites involved in reward, affect, social cog- 
nition, and emotion regulation, particularly the midbrain, stria- 
tum, amygdala, hippocampus, anterior cingulate cortex (ACC), 
and orbitofrontal cortex (OFC) (Olazabal and Young, 2006; 
Campbell, 2008). It has been hypothesized that the upregu- 
lation of OT receptors during late pregnancy and parturition 
leads directly or indirectly, via dopaminergic reward circuits, 
to reinforcing interactions with offspring, reduced avoidance of 
offspring, adaptive maternal behavior, and increased protective 
maternal aggression against intruders in small-brain mammals 
(Numan and Woodside, 2010). 

Given this literature we focus our attention on a gene with 
conceptual linkages to caregiving behaviors and parenting. The 
human OXTR gene is located on chromosome 3p25, contain- 
ing four exons and three introns. An SNP in the third intron of 
OXTR, rs53576 (G/A), has been identified as a candidate gene 
underlying social behavior in humans (Inoue et al., 1994). Only 
a small literature exists on associations of OXTR gene polymor- 
phisms with human behavior. Individuals with the TT genotype 
of the rsl042778 SNP in the 3' untranslated region of OXTR 
allocated more resources to an unseen opponent in a competi- 
tive game (Israel et al., 2009). Individuals with the GG genotype 
of the intronic SNP rs53576 in OXTR exhibited higher disposi- 
tional empathy and less stress reactivity (Rodrigues et al, 2009), 
whereas the AA genotype was associated with slightly higher 
depression and lower self-esteem scores (Saphire-Bernstein et al., 
201 1). In another study, however, the AA genotype of rs53576 was 



associated with lower positive affect in males only and not related 
to negative affect or loneliness (Lucht et al., 2009). 

Tost et al. (2010) used a partial genetic imaging paradigm to 
conduct tests of associations between genetic polymorphisms and 
both neural and behavioral phenotypes. The rs53576 GG geno- 
type was associated with greater dependence on the approval 
of others and greater amygdala activation and lower amygdala- 
hypothalamus coupling in functional magnetic imaging (fMRI) 
in a face-matching task (Tost et al., 2010). The GG genotype 
was associated with smaller hypothalamic volume, which in turn 
was correlated with dependence on the approval of others, but 
predominantly in males. No tests of associations between brain 
regions activated in fMRI and dependence on social reward were 
reported, however. 

For purposes of the present study, the existing literature has 
limitations. Most studies used college students (Israel et al., 2009; 
Rodrigues et al., 2009; Saphire-Bernstein et al, 2011) or volun- 
teers screened to exclude mental disorders (Tost et al, 2010), 
minimizing both women with children and participants with 
impaired psychological functioning. Interestingly, the findings of 
the one study of a representative sample were inconsistent with 
the studies of less representative samples (Lucht et al, 2009). 
Additionally, there is limited evidence that the behavioral phe- 
notypes used in these studies are correlated with variations in 
maternal parenting, except for depressive symptoms and neg- 
ative affect, which have only been limitedly linked to rs53576 
(Saphire-Bernstein et al., 2011). 

A few studies of the relevance of OT to human maternal par- 
enting have been published recently, however. One found that 
plasma OT during pregnancy and the postpartum was positively 
correlated with adaptive maternal parenting (Levine et al., 2007). 
In addition, intranasal OT decreased amygdala activation and 
increased activation of the insula and inferior frontal gyrus when 
women heard infant cries (Riem et al., 201 la). 

Directly relevant to the present study, a study of Caucasian 
mothers of 2-year old children with behavior problems found that 
mothers with the OXTR rs53576 GG genotype displayed more 
sensitive parenting (Bakermans-Kranenburg and van Ijzendoorn, 
2008). In a study with more complex findings, however, nulli- 
parous adult females with the rs53576 GG genotype had greater 
heart rate responses to infant cries, but only among women with 
low depression scores, with the opposite finding for depressive 
women (Riem et al., 201 lb). 

To further lay a foundation for large-scale studies of the neu- 
robiology of dysfunctional maternal parenting, we conducted 
a study of mothers selected based on their scores on positive 
and negative parenting of their 4-6 year old children. Using a 
full genetic imaging paradigm, we tested hypotheses that OXTR 
gene polymorphisms are associated with variations in both neu- 
ral responses to child stimuli and maternal parenting, and that 
neural responses in the same regions are correlated with mater- 
nal parenting. fMRI studies support such genetic imaging stud- 
ies of human maternal parenting by showing that infant faces 
and cries activate maternal brain regions shown to contain OT 
receptors in rodents (i.e., amygdala, hippocampus, and striatum) 
(Lee et al., 2009; Stoop, 2012), and midbrain and basal ganglia 
regions involved in reward (Lorberbaum et al., 2002; Bartels and 
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Zeki, 2004; Swain et al., 2007). In more experienced mothers 
of older children, amygdala and insula activations are not typi- 
cally observed, however, but regions involved in social cognition, 
such as the OFC, medial prefrontal cortex (MPFC), and supe- 
rior temporal sulcus (STS), are activated (Leibenluft et al., 2004; 
Swain et al, 2004; Swain, 2011). No fMRI studies have exam- 
ined associations between OXTR gene variants and variations in 
both maternal parenting and maternal neural responses to child 
stimuli in a genetic imaging paradigm. 

METHODS 
PARTICIPANTS 

Seventeen female and 104 males 4-6 year olds were recruited 
from a child psychiatry clinic in Chicago during 1994-1995 for a 
longitudinal study of children with attention-deficit/hyperactivity 
disorder (ADHD) and matched controls (Lahey et al., 2005; 
Chronis-Tuscano et al., 2010; Lee et al, 2010). Participants tak- 
ing short-acting stimulant medications were included, but all 
assessments were conducted when the child was off medication. 
Previous reports from this study revealed that the measures of 
positive and negative maternal parenting used predicted poor 
outcomes of children with ADHD through adolescence, control- 
ling for multiple factors (Chronis et al, 2007; Lahey et al., 2011). 
Thus, the genetic and neural correlates of measures of maternal 
parenting already known to predict long-term child outcomes 
can be studied in this sample. In 2009-2010, 40 mothers were 
selected for the present study by recruiting from the extremes of 
positive and negative parenting scores to maximize variation in 
parenting. Participants' written consent was obtained. All partic- 
ipants were paid for their participation. The study was approved 
by the University of Chicago Institutional Review Board and con- 
ducted in accordance with the Declaration of Helsinki. All of the 
observed children had written informed consent from their bio- 
logical mothers, all of whom were custodial. Characteristics of 
participants are in Supplemental Table SI. 

MOTHER-CHILD INTERACTION TASK 

Mother-child interactions were videotaped in a room equipped 
with furniture, toys and other objects, and a television showing 
cartoons (Chronis et al., 2007). Mothers were invited to play freely 
with their children. After lOmin, the interviewer returned and 
scattered clothes, papers, and empty containers around the room. 
The mother was given an Etch-a- Sketch, worksheets, a magazine, 
a pencil, and written instructions to complete tasks with her child 
in order over 15min: (1) return toys to shelves, (2) put clothes 
in the box, (3) put paper and containers in the wastebasket, (4) 
count geometric shapes, (5) copy a set of geometric designs on 
paper, (6) dust the table with a cloth, (7) draw a diagonal line on 
the Etch-a-Sketch, and (8) choose one toy and play quietly while 
the mother reads and takes a telephone call. The first 13min of 
the task situation were coded because some dyads completed all 
tasks early. 

Interactions were coded using the Dyadic Parent-Child 
Interaction Coding System (DPICS) (Robinson and Eyberg, 
1981). Two reliable and valid measures of parenting were used 
by averaging standardized scores across the structured and play 
situations: Positive parenting (praise, positive affect, and physical 



positive) and negative parenting (negative commands, critical 
statements, and physical negatives) (Robinson and Eyberg, 1981; 
Webster- Stratton, 1998). These measures predict long-term child 
outcomes in this sample (Chronis et al, 2007; Lahey et al., 2011). 
Thirty percent of the videotapes were coded by a reliability coder. 
Inter-observer agreement was 0.95 for positive parenting, 0.90 for 
negative parenting, and 0.92 for total child disruptive behavior. 

NEUR0IMAGING 

Photographs of the child at 4-6 years from the videotapes and 
photos of unrelated demographically matched children were used 
in fMRI as stimuli. Mothers also viewed a series of dynamic 
visual stimuli, each consisting of three 600 x 480 pixel color 
photographs presented successively for 1000, 200, and 1000 ms, 
respectively to imply motion. Forty-eight dynamic stimuli por- 
trayed misbehaviors (e.g., a child intentionally kicking a female 
adult on the leg), and 48 portrayed parallel neutral behaviors (e.g., 
a child standing next to a female adult) without showing faces. 
The stimuli were matched on numbers of people and objects and 
varied in skin color. 

Stimuli were presented with E-prime 1.2 (PsychologySoftware 
Tools, 201 1) by back-projection in a block design. Dynamic child 
stimuli were presented in 16 19.8 s blocks and a fixation cross was 
presented in 16 18 s baseline blocks. Stimuli were blocked by type 
(provocative/neutral child behavior), each consisting of 6 stimuli 
(2200 ms each) with six 1100 ms inter-stimulus intervals, during 
which a black fixation cross was presented against a gray back- 
ground. Participants were shown the stimuli in 2 sessions (8 active 
blocks per session). In one session, immediately preceding every 
stimulus block, mothers were shown a photograph of their own 
child for 6 s and instructed to "imagine this is your child" in the 
actions that followed. In the other session, mothers were shown 
a photograph of a demographically matched child and instructed 
to "imagine this is not your child" in the actions that followed. 
Participants were instructed to watch the stimuli carefully and no 
overt response was required. Session order was counterbalanced 
across participants. 

MRI was performed on a 3-T Philips Achieva Quasar scan- 
ner. The fMRI pulse sequence parameters include time repeti- 
tion/time echo (TR/TE) 2000/25, flip angle = 77, 32 contiguous 
slices with 4 mm thickness, slice gap 0.5 mm, 224 x 224 mm 2 
field of view (FOV), approximately 64 x 64 matrix. High res- 
olution structural images were acquired in the sagittal plane 
using a Tl-weighted 3D Turbo Field Echo (TFE/MP-RAGE) 
anatomical scan with the following parameters: ri? = 8.1ms, 
TE = 3.7 ms, FOV = 224 x 224 x 160 mm 3 , isotropic voxel size 
lxlxl mm 3 , matrix size 224 x 224. During anatomical scans 
after the stimulation paradigms, participants watched a movie 
about tropical beaches. 

GEN0TYPING 

DNA from mothers was isolated from saliva in Oragene kits 
(DNAGenotek, 2010) and checked for quality by OD ratio of 
260/280 and concentration. Genotyping rs53576 and rsl042778 
was performed using TaqMan pre-designed SNP genotyping 
assays (AppliedBiosystems, 2010). PCR was carried out in a total 
volume of 3 |xl containing 10 ng genomic DNA, 1.5 (Jtl 2 x TaqMan 
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universal PCR master mix (AppliedBiosystems, 2010), 0.075 ill 
40 x SNP genotyping assay. After 95° C lOmin, 40 cycles consist- 
ing of 15 s at 92° C and 1 min at 60° C annealing temperature were 
performed. After PCR amplification, an endpoint plate read using 
a 7900 Real-Time PCR System was performed. The Sequence 
Detection System (AppliedBiosystems, 2010) was used to call 
genotypes. 

STATISTICAL ANALYSES 
fMRI analysis 

Image processing was carried out with SPM8 in MATLAB 7.0 
(Marsh et al., 2008). Preprocessing included correction for head 
motion, normalization to the SPM8 echo-planar imaging tem- 
plate, and smoothing using a 6-mm full-width half-maximum 
isotropic Gaussian kernel. Images were realigned and normal- 
ized using standard SPM procedures. All participants had less 
than 0.5 voxels of in-plane motion throughout scanning. A 2-level 
approach for block-design fMRI data was adopted using SPM8. 
A voxel-by-voxel multiple regression analysis of expected signal 
changes for the child photographs and the 2 block categories, con- 
structed using the SPM8 hemodynamic response function, was 
applied to preprocessed images. Individual subject data were ana- 
lyzed using a fixed-effects model. Group data were analyzed using 
a random-effects model. Condition effects at the subject level 
were modeled by box-car regressors representing type of child 
photograph and the two block types. 

The fMRI contrasts that were examined across participants 
were: ( 1 ) Own vs. other child image preceding each of the action 
blocks and (2) provocative vs. neutral child behavior in the 
own-child session. Activations were overlaid on a representative 
high-resolution structural Tl -weighted image from one subject 
from the SPM8 canonical image set, coregistered to Montreal 
Neurological Institute (MNI) space. 

To additionally test for activations previously reported in 
maternal parenting studies, we selected five regions of interest 
(ROIs) identified in previous published work on the neuroscience 
of maternal parenting (Swain, 2011) and about which we had a 
priori hypotheses, using a small volume correction (SVC) for fam- 
ily wise error (FWE) at P < 0.05. The following bilateral regions 
were selected: amygdala, hippocampus, anterior insula, ACC, 
and OFC. Analyses were implemented in SPM8 based on cor- 
rections for multiple comparisons confined to an ROI (Worsley 
et al., 1996). Data extraction for 5-mm spherical ROIs was per- 
formed using the rfxplot toolbox (Glascher, 2009) in SPM8. 
Coordinates were based on results from the whole-brain anal- 
yses and neuroanatomical atlases. The small volumes consisted 
of a 5 mm sphere centered at the most significant voxel of the 
clusters activated at P < 0.001 uncorrected in the whole brain 
analysis. 

Tests of association in genetic imaging paradigm 

Three sequenced and prioritized sets of generalized linear models 
(Nelder and Wedderburn, 1972) were conducted as illustrated in 
Supplemental Figure SI: 

(1) Additive terms were coded —1, 0, and 1 to jointly test lin- 
ear differences among 0, 1, or 2 minor A alleles for rs53576 



and G alleles for rsl042778. Non-additive terms, coded —1, 
2, and —1, jointly captured any mode of non-additive trans- 
mission for each SNP. Additive and non-additive terms were 
run in the same model. In two separate analyses, positive 
and negative maternal parenting were regressed on the four 
OXTR terms (additive and non-additive terms for rs53576 
and rsl042778), simultaneously controlling child's sex, race- 
ethnicity, birth order, age in wave 1, mother's age at scanning, 
delivery (Caesarian or vaginal), and total parity. To adjust for 
child effects on the mother's parenting (Bell and Chapman, 
1986), the child's diagnosis of ADHD and disruptive behavior 
during the mother-child interaction also were simultaneous 
covariates (Lahey et al., 2011). 

(2) For each parenting measure for which the tests of genetic 
associations were significant atp < 0.05 in step 1, generalized 
linear models regressed parenting on each of the five bilateral 
a priori selected ROIs, with the same covariates plus time of 
fMRI scanning. 

(3) For each of the five a priori selected ROIs that was signifi- 
cantly associated with parenting in step 2, regression models 
with the same covariates as step 2 tested associations between 
the SNP and the ROI. 

RESULTS 

WHOLE-BRAIN ANALYSES 

In the whole-brain analyses, the own > other child contrast was 
associated with hemodynamic increases in regions subserving 
motivation, reward, and emotion regulation processing (includ- 
ing the midbrain, dorsal putamen, thalamus, anterior cingulate, 
and prefrontal cortices) and previously been found to be acti- 
vated by child stimuli (Swain, 201 1). The inappropriate > neutral 
behavior contrast activated similar regions, with additional acti- 
vations in areas involved in social cognition such as the posterior 
STS and deactivations in somatosensory cortex and bilateral 
posterior insula (Figures 1, 2 and Supplemental Tables S2, S3). 

ASSOCIATIONS BETWEEN SNPs AND PARENTING 

Allele frequencies are shown in Supplemental Table S4. There 
were no significant ancestry group differences in allele frequen- 
cies for rs53576, but the T allele of rsl042778 was more common 
in African Americans, requiring ancestry to be covaried. There 
were no significant associations between OXTR SNPs and neg- 
ative parenting, but rs53576 and rs 1042778 each accounted for 
unique variance in positive maternal parenting, although only the 
additive association between rs53576 alleles and positive parent- 
ing was significant atp < 0.05/8 = 0.006 corrected for multiple 
testing (Table 1). A subsequent regression analysis revealed sig- 
nificant interactions between ancestry and rs53576 (non-additive 
term: x 2 = 6.05, P < 0.02) and rsl042778 (additive term: \ 2 = 
10.11, P < 0.002). Although cell sizes were too small to interpret 
these interactions, they raise the possibility that OXTR polymor- 
phisms may be associated with parenting differently in mothers 
of African and European descent. 

ASSOCIATIONS BETWEEN ROIs AND PARENTING 

In separate analyses for each ROI, we observed associations 
between observed positive parenting and activation in bilateral 
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FIGURE 1 | Maternal hemodynamic brain activations in response to photographs of her own vs. an unrelated child (p - 
Supplemental Table S2 for corrected activations. 



0.005, uncorrected). See 



Hippocampus 



Somatosensory 
Cortex 







FIGURE 2 | Maternal hemodynamic brain activations in response to dynamic images of a young person behaving inappropriately vs. appropriately 
(p < 0.005, uncorrected). See Supplemental Table S3 for corrected activations. 



OFC and left ACC when mothers viewed their own vs. another 
child, at P = 0.05 levels (see Table 2). Mothers who scored higher 
on observed positive parenting behaviors showed an increase 
in hemodynamic response in OFC and ACC when they viewed 
pictures of their own child, compared to pictures of another 
child. 

Two associations between ROIs activated in the inappropri- 
ate > appropriate child behavior contrast and positive parenting 
were significant at P = 0.05 levels. Associations were observed 
in the right ACC and right hippocampus (see Table 2). Because 
we focused our analyses on specific regions, hypothesized a pri- 
ori, ROIs above aP = 0.05 uncorrected threshold were included 
in the next step. It should be noted, however, that none of the 
analyses were significant at a = 0.05/20 = 0.0025 corrected for 
multiple testing. 

ASSOCIATIONS BETWEEN SNPs AND ROIs 

As shown in Figure 3 and Table 3 tests of additive associations 
between rs53576 and the ROIs that were both significantly 
activated by child stimuli and were significantly associated with 
positive parenting, were significant at p < 0.05 uncorrected 



levels. Specifically, these were in bilateral OFC, left ACC and 
right hippocampus. Further, the association between left 
ACC activation and rs53576 (additive) in the own > other 
child contrast and between right hippocampus activation and 
rs53576 (additive) in the inappropriate > appropriate behavior 
contrast remained significant after Bonferroni correction 
(a = 0.05/20 = 0.0025). 

In addition, one test of independent associations between 
rsl042778 (non-additive) and these ROIs were significant atp < 
0.05 levels (uncorrected). This was in the right hippocampus. 
None were significant at a = 0.05/20 = 0.0025 levels. 

DISCUSSION 

Consistent with previous studies of human mothers (Lorberbaum 
et al., 2002; Nitschke et al., 2004; Swain et al, 2007, 2008), 
two kinds of child stimuli activated neural regions impli- 
cated in face processing, reward processing, social cogni- 
tion, and executive control of emotion (Figures 1, 2 and 
Supplemental Tables S2, S3). More importantly, individual dif- 
ferences in maternal neural responses in key ROIs in both 
paradigms were associated with variations in maternal parenting. 
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Notably, the present findings are internally consistent. Not only 
was the rs53576 A allele associated with higher levels of positive 
parenting, it also was associated with greater activation in the 
own > other child contrast in two a priori hypothesized ROIs 
that were associated with positive parenting (Figure 3). In par- 
ticular, the rs53576 A allele was associated with greater activation 
in ROIs involved in (1) orienting toward, monitoring, and evalu- 
ating child cues (Parsons et al., 2013) and other emotional stimuli 
(bilateral OFC) (Rolls, 2000; Nitschke et al, 2004; Volkow et al, 
2011) and (2) assessing the salience of emotional information 



and regulating emotional responses (left ACC) (Bush et al., 2000; 
Kerns et al., 2004). The present findings are consistent with previ- 
ous findings that OFC activation is correlated with positive mood 
related to maternal behavior (Nitschke et al., 2004). Furthermore, 



Table 2 | Tests of associations between observed maternal positive 
parenting and maternal BOLD activations by the two kids of stimuli 
in brain regions of interest using small-volume correction (N = 34). 



Cortical regions of interest 



Positive Parenting 



Table 1 | Simultaneous associations of each term for additive and 
non-additive associations (each 1 DF) between OXTR SNPs with 
positive and negative maternal parenting (in separate analyses), 
adjusting for all covariates; N = 35. 



POSITIVE PARENTING 

rs53576 (additive) 
rs53576 (non-additive) 
rs1042778 (additive) 
rs1042778 (non-additive) 
NEGATIVE PARENTING 
rs53576 (additive) 
rs53576 (non-additive) 
rs1042778 (additive) 
rs1042778 (non-additive) 



15.43 
0.02 
2.93 
5.37 



0.09 
0.01 
1.00 
0.01 



Wald 



P = 



<0.0001 
0.8814 
0.0872 
0.0205 



0.7647 
0.9339 
0.3178 
0.9305 



Simultaneous covariates: Child's age in years, child's sex, African American 
ancestry, child's ADHD diagnosis, child's disruptive behavior during mother-child 
interaction task, first born status, mother's number of live births, prematurity in 
weeks, Caesarian delivery, and mother's age at scanning. 



OWN > UNRELATED CHILD 
Orbital frontal cortex 

Left 
Right 



0.25 
0.26 



4.56* 
4.60* 



Anterior cingulate cortex 

Left 0.43 5.64* 

Right 0.12 0.44 

PROVOCATIVE > NEUTRAL BEHAVIOR WHEN VIEWING OWN CHILD 

CONTRAST 

Anterior cingulate cortex 

Left -0.15 0.10 

Right -0.97 5.87* 



Hippocampus 

Left 
Right 



0.18 
1.42 



0.08 
3.85* 



*P < 0.05. Simultaneous regression analyses controlling for the child's sex and 
age at the time of the mother-child interaction task, the child's race-ethnicity 
the child's diagnosis (ADHD or healthy control), the child's level of disruptive 
behavior during the task, premature birth, caesarian delivery, parity, birth order, 
time of day of the mother's scan, and mother's age at the time of the scan. 
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FIGURE 3 | Scatter plots with group means residualized on all covariates 
for hemodynamic responses in the own/other child (O) contrast and the 
inappropriate/appropriate behavior contrast (I) in brain regions of 
interest that were found to be significantly related to oxytocin and 



vasopressin genotypes, and for observed maternal parenting, presented 
by oxytocin receptor gene (OXTR) rs53576 genotypes. To facilitate 
interpretation, all variables on y axes are expressed as sample z-scores 
(mean = 0; standard deviation = 1). 
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Table 3 | Results of simultaneous regression analyses testing 
associations of each additive and non-additive genetic term (each 
DF — 1) for OXTR rs53576 and rs1042778 with each neural region of 
interest that was both significantly activated in the fMRI paradigm 
and found to be significantly associated with observed positive 
maternal parenting with covariates, conducted separately for each 
neural region of interest (/V = 30). 

Wald x 2 P < 



OWN CHILD > OTHER CHILD fMRI PROBE PARADIGM 

Left orbital frontal cortex 

rs53576 (additive) 3.90 

rs1042778 (non-additive) 0.06 

Right orbital frontal cortex 

rs53576 (additive) 4.67 

rs1042778 (non-additive) 0.00 

Left anterior cingulate cortex 

rs53576 (additive) 29.07 

rs1042778 (non-additive) 0.38 



INAPPROPRIATE > NEUTRAL CHILD BEHAVIOR fMRI PROBE 
PARADIGM 

Right anterior cingulate cortex 

rs53576 (additive) 

rs1042778 (non-additive) 
Right hippocampus 

rs53576 (additive) 

rs1042778 (non-additive) 



0.10 
0.26 



17.42 
5.69 



0.0482 

0.8080 

0.0306 

0.9806 

<0.0001 

0.5398 



0.7552 
0.6113 

0.0001 
0.0171 



Simultaneous covariates in each model: child's age in years, child's sex, African 
American ancestry, child's ADHD diagnosis, child's disruptive behavior during 
mother-child interaction task, first born status, mother's number of live births, 
prematurity in weeks, Caesarian delivery, time of day of scanning, and mother's 
age at scanning. 



when mothers viewed inappropriate behavior when instructed to 
imagine it was their own child, the rs53576 A allele was associated 
with greater activation of the hippocampus. This is interesting 
given that activation of OT receptors in the hippocampus in rats 
inhibits behavioral responses to stressors (Cohen et al., 2010). 
Because child misbehavior often elicits harsh maternal behavior 
(Lahey et al, 1984, 2011), this may reflect greater inhibition of 
such responses in more positive mothers. 

Mothers with T alleles of OXTR rsl042778 exhibited both 
greater positive parenting and greater activation of the right hip- 
pocampus in the own > other child contrast, but otherwise 
there was less clarity and internal consistency in the findings for 
rsl042778. Furthermore, the nature of the non-additive associa- 
tions of rsl042778 with the hippocampal ROI, and its relation to 
positive parenting, is unclear. 

Nonetheless, it is essential to note that the present study 
replicated the association between genotype and positive parent- 
ing found in the only previous study of rs53576 and parenting 
(Bakermans-Kranenburg and van Ijzendoorn, 2008) in reverse. In 
the previous study, mothers with GG genotype of rs53576 dis- 
played more sensitive parenting with 2 year olds than mothers 
with the AA or AG genotype. Because such "allele flip" replica- 
tions are commonly reported in the genetics literature, Clarke 



and Cardon (2010) delineated the conditions under which the 
results of both studies could be valid. This is possible when effect 
sizes are large, the samples differ in ancestry composition, and 
the target polymorphism interacts with unmeasured polymor- 
phisms that vary in allele frequencies among ancestry groups 
(Clarke and Cardon, 2010). Because the first two conditions are 
met in this case, interactions of both SNPs with other polymor- 
phisms that differ in allele frequency in persons of African and 
European descent could result in the differences in findings for 
rs53576 in the two studies. Thus, the present findings could be a 
valid "allele flip" replication of the previous studies (Bakermans- 
Kranenburg and van Ijzendoorn, 2008). No conclusions can be 
reached, of course, until the hypothesis that the direction of asso- 
ciations between genotypes of rs53576 and parenting differ by 
ancestry group are tested in the future. 

One potential difference between these two studies is that 
mothers in the present study were rated on positive parenting 
behaviors, whereas the study by Bakermans-Kranenburg and van 
Ijzendoorn assessed sensitive parenting, which may be separable 
constructs. A consistently positive mother, for instance, may not 
always be sensitive to her child's needs, while a sensitive mother 
may not always engage in positive parenting behaviors. 

It is important to understand the present findings in the 
context of several issues and limitations: 

LIMITED MARKERS 

Only two SNPs were used to characterize variation in the 
OXTR gene. Thus, although the present findings support further 
research, future studies should genotype all polymorphisms in 
this and related genes to fully understand their impact on neural 
responses to child stimuli and variations in maternal parenting. 
This would likely require a larger sample size. 

CONSERVATIVE STATISTICAL STRATEGY 

The statistical strategy illustrated in Supplemental Figure SI lim- 
ited the statistical analyses of hemodynamic responses to child 
stimuli to only those regions that were both significantly activated 
by child stimuli and were significantly associated with positive 
parenting, and only included SNPs that were found to be sig- 
nificantly associated with positive parenting. This was done to 
focus on the feasibility of testing hypotheses regarding neural 
pathways from variations in these SNPs to variations in mater- 
nal parenting using fMRI-based genetic imaging methods. Thus, 
for example, tests were not conducted to determine if OXTR gene 
polymorphisms were associated with hemodynamic responses to 
child stimuli in ROIs that were not significantly associated with 
maternal parenting. Thus, these analyses could underestimate the 
number of ROIs whose functioning is associated with these SNPs. 

SAMPLE SIZE 

In principle, the size of a sample is unrelated to the likelihood of 
rejecting the null hypothesis when there is no association in the 
population (type I error). This is because the alpha level (e.g., p < 
0.05) for each individual statistical test is set in advance and the 
statistical test is automatically penalized for smaller sample sizes. 
Therefore, the probability that the null hypothesis was rejected 
incorrectly is not a major issue in the present study. Nonetheless, 
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smaller samples impose three serious limits on statistical testing. 
First, because the statistical tests used in this study are based on 
estimations, a degree of error can be expected in smaller sam- 
ples. Second, the probability of failing to reject the null hypothesis 
when there is an association in the population (type II error) is 
jointly related to the sample size and the magnitude of the effect 
size. Only associations with large effect sizes can be detected in 
small samples. Therefore, it is possible that the present study 
failed to detect associations of smaller but still important effect 
sizes. Third, confidence intervals around effect sizes are larger in 
smaller samples. This makes it more difficult to determine the 
sample size needed in future studies. This is because power analy- 
ses must be based on the lower limit of the confidence interval for 
the effect size and smaller sample sizes provide lower estimates of 
this value. 

SAMPLE REPRESENTATIVENESS AND HETEROGENEITY 

A potentially more serious limitation is that the present sample is 
not representative, partly because smaller samples are inherently 
less likely to represent the population. In addition, this sample 
was selected from a single clinic rather than from a representa- 
tive sampling frame. This means that estimates of the magnitude 
of associations are unlikely to accurately mirror the population. 
Furthermore, the present sample also is heterogeneous in terms of 
race-ethnicity and in the different methods used to sample cases 
between ADHD and controls. Although a range of demographic 
variables were controlled in the present analyses, it is possible that 
the findings are influenced to an unknown degree by artifacts of 
population stratification of alleles. Future studies must address 
these issues by using samples that are either more homogeneous 
or more representative and large enough to adjust for population 
stratification. 

TIME BETWEEN OBSERVATION OF MOTHER-CHILD INTERACTION AND 
SCANNING 

In this initial study, we used an available sample of mothers 
on whom mother-child interactions had been observed when 
the children were 4-6 years old. This allowed us to select the 
sample for scanning and genotyping by oversampling on the 
extremes of maternal parenting in a way that ensured a broad 
distribution of maternal parenting scores to increase statistical 
power for detecting associations with parenting. For this particu- 
lar sample, this necessitated scanning the mothers as they viewed 
stimuli approximately 15 years later. The implications of this time 
frame must be considered separately for each set of statistical 
tests. 

The 15 -year gap has no implications for the tests of associ- 
ation between genotypes and hemodynamic responses of ROIs, 
as genotyping and scanning were conducted contemporaneously. 
The gap similarly is not likely to be an issue for the tests of associa- 
tions between genotypes and parenting as DNA sequences do not 
change over time. Nonetheless, it is possible that our tests of asso- 
ciation between parenting and genotypes are conservative. This 
is because epigenetic changes that affected expression of those 
genotypes could have occurred during the intervening 15 years. 
The gap between the mother-child observations and scanning 
could have weakened or even distorted these tests of associations 



between ROI responses and parenting. For instance, the differ- 
ence in ROI activation between genotypes may potentially reflect 
an effect of genotype on how mothers remember their child or 
their child's behavior, or even the status of the current relation- 
ship with the child. Relatedly, given this 15 -year gap, the mother's 
current relationship with her child might impact any correlations 
between positive parenting and hemodynamic activity in ROIs, 
particularly as the parenting measures were previously found to 
predict long-term outcomes for this sample (Chronis et al., 2007; 
Lahey et al., 2011). Nonetheless, the time gap had the effect of 
focusing the present study on only highly enduring aspects of 
neurobiological functioning that are related to variations in early 
maternal parenting. It is also plausible that the perceived effects 
of genotype on behavior and/or brain response might be due to 
the presence of the same or similar genotype in either the previ- 
ous generation (possibly affecting the mother's own upbringing) 
or in the child (possibly affecting the behavior of the child and, 
therefore, the associations and responses the mother has toward 
them). 

The ultimate goal of this line of research is to trace environ- 
mentally moderated neural pathways from variants in genes to 
dysfunctional maternal parenting to improve prevention. Such 
research is essential to understanding why some mothers parent 
adaptively whereas others are neglectful or abusive. As such, this 
study was designed to provide a foundation for future research, 
rather than to conduct strong tests of hypotheses based on previ- 
ous studies. 

To our knowledge, the present study is the first to provide 
linked evidence on such putative pathways from polymorphisms 
in OXTR to neural functioning to parenting. These findings sup- 
port future tests of the hypotheses that variations in genes coding 
OXTR are related to variations in human maternal parenting. 
Although the present findings are internally consistent, they are 
not consistent with a previous study and there are limitations 
to this study regarding the sample and the time between the 
observed mother-child interactions and scanning. 
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00021 /abstract 

Supplemental Figure SI | Flowchart of the sequenced statistical analyses. 

Supplemental Table S1 | Demographic and birth characteristics of three 
subgroups of a total 40 mothers and their children participating in tests of 
associations between different variables. 

Supplemental Table S2 | Areas of significant hemodynamic activation 
when viewing a picture of own child > unrelated child. 
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Supplemental Table S3 | Areas of significant activation and deactivation 
when viewing a person behaving inappropriately > behaving neutrally in 
the own child condition. 

Supplemental Table S4 | Allele frequencies of OXTR SNPs by ancestry 
groups. 
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